AD-A264  733 


NSWCDD/TR-92/243 


PHALANX  CIWS  CONTROL  SYSTEM  STABILITY,  AIM 
BIAS  COMPENSATION,  AND  NOISE  SENSITIVITY 


BY  DEMETRIOS  SERAKOS 
COMBAT  SYSTEMS  DEPARTMENT 


MAY  t992 


Approved  for  public  release,  distribution  is  unlimited. 


NAVAL  SURFACE  WARFARE  CENTER 

DAHIGREN  DIVISION 

Dahlgr«n.  Virginia  22448-5000 


93-11286 


NSWCDD/TR-92/243 


PHALANX  CIWS  CONTROL  SYSTEM  STABILITY.  AIM 
BIAS  COMPENSATION,  AND  NOISE  SENSITIVITY 


BY  DEMETRIOS  SERAKOS 
COMBAT  SYSTEMS  DEPARTMENT 


MAY  1992 


Approved  for  public  release;  distribution  is  unlimited 


NAVAL  SURFACE  WARFARE  CENTER 
DAHLGREN  DIVISION 
Oahlgren,  Virginia  22448-5000 


NSWCDD/TR-92/243 


FOREWORD 


This  report  presents  a  stability  analysis  of  an  aiming  control  system  that  is  similar  to  that  used  in 
the  Block  I  PHALANX  Oose-In  Weapon  System  (CIWS).  Besides  stability,  other  design  goals  of  the 
control  system  are  aim  bias  compensation  and  noise  sensitivity.  This  report  presents  a  theoretical 
analysis  of  these. disparate  goals  and  design  tradeoffs  are  also  presented.  It  is  a  continuation  of  previ¬ 
ous  work  reported  in  Bailey,  E.P.  and  Price  E.L.,  An  Analysis  of  Gun  Aim  Bias  Compensation  for  the 
PHALANX  CIWS,  FMC  Corporation,  King  George.  Virginia.  May  1987,  which  largely  was  an  analysis 
of  improved  performance  over  a  previous  control  system.  Control  theoretic  issues  such  as  stability, 
aim  bias  compensation  and  noise  sensitivity  were  not  considered  in  that  work. 

The  author  wishes  to  thank  E.L.  Price  and  M.H.  Pee  of  FMC  Corporation  for  introducing  this 
subject  and  for  their  technical  advice. 

This  report  has  been  reviewed  by  A.  Garza.  Head,  System  Engineering  Branch  and  R.E.  Lutman, 
Head,  AEGIS  Ship  Combat  Systems  Division. 
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ABSTRACT 


An  aiming  control  system,  which  is  similar  to  that  in  the  Block  I  PHALANX  Close  In  Weapon 
System  is  considered  in  this  report.  An  important  feature  of  this  control  system  is  that  it  compensates 
for  any  gun  aim  bias.  An  aim  bias  may  be  caused  by  variations  in  the  gun.  ammunition,  or  environ¬ 
mental  conditions.  Design  issues  considered  are  stability,  aim  bias  compensation,  and  sensitivity 
feedback  noise.  These  are  disparate  design  goals.  Design  tradeoffs  ’hat  quantify  this  disparity 
presented.  An  example  is  given  that  illustrates  hew  the  analysis  developed  in  this  report  might 
used  in  a  design  situation. 
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INTRODUCTION 


In  this  report  an  aiming  control  system  that  is  similar  to  that  in  the  Block  I  PHALANX  Qose-In 
Weapons  System  (CIWS)  is  considered.  Among  the  important  issues  considered  in  designing  this 
control  system  are  stability,  target  tracking,  distuihance  rejection  and  sensitivity  to  noise.  As  it  turns 
out,  designing  a  system  with  better  target  tracking  makes  it  less  stable.  Also,  better  target  tracking 
makes  the  system  more  sensitive  to  noise.  A  system  less  susceptible  to  noise  will  be  more  susc  jf)i-ble 
to  disturbances.  Previously,  no  analytical  work  has  been  done  quantifying  these  design  considerations. 
The  purpose  of  this  report  is  to  present  an  analytical  picture  to  these  design  considerations  that  may  be 
used  to  sharpen  a  current  design  or  speed  the  development  of  a  future  design. 

A  broad  picture  of  the  problem  is  as  follows:  The  desired  gun  aiming  angle  is  fed  into  the  con¬ 
trol  system.  The  aiming  angle  is  treated  as  a  reference  input  -  the  manner  in  which  it  is  obtained  is 
not  considered  in  this  report.  The  gun  tracks  this  reference  input.  A  proportional  plus  integral  con¬ 
troller  is  used  in  an  inner  feedback  loop  to  accomplish  this.  One  problem  is  that  there  may  be  an  aim 
bias.  Among  the  factors  contributing  to  an  aim  bias  are:  misalignment  of  the  gun  aiming  mechanism, 
the  barrels  heating  due  to  the  gun  being  fired,  barrel  wear,  variations  in  the  propellant  used,  wind, 
temperature,  and  barrel  warping  due  to  sun  load.  The  aim  bias  is  modeled  as  a  disturbance  input. 
The  control  system  features  a  second  outer  feedback  loop  to  null  out  this  bias.  The  bullet  stream 
angle  at  the  target  is  measured.  Ideally,  this  stream  should  match  the  measured  gun  angle  at  the  time 
those  bullets  were  fired.  If  the  two  angles  do  not  match,  the  error  angle  is  fed  back  around  the  outer 
loop  to  compensate.  The  bullet  stream  angle  measurement  is  noisy,  so  a  noi,se  filter  is  included  in  the 
outer  feedback  loop. 

In  the  PHALANX  CIWS  Aiming  Model  section,  a  model  for  this  system  is  defined  and  dis¬ 
cussed.  The  transfer  functions  that  are  needed  arc  computed.  It  will  be  seen  that  this  is  a  linear  sys- 
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tern  with  a  time-varying  time  delay.  Time  delays  are  destabilizing.  The  time  delay  arises  because  the 
feedback  to  null  out  the  aim  bias  involves  past  reference  angles.  The  stability  of  the  system  is  then 
considered.  The  first  consideration  is  the  stability  of  the  inner  loop.  This  part  of  the  system  is  a 
linear  system  (with  no  time  delays).  The  general,  but  conservative,  Small-Gain  Theorem',  is  used  to 
give  conditions  for  the  closed-loop  stability  of  the  overall  system.  The  Small-Gain  Theorem  applies  to 
systems  with  time-varying  time  delays.  Roughly  speaking,  the  Small-Gain  Theorem  states  that  a  feed¬ 
back  system  is  stable  if  the  gain  of  the  transfer  functions  around  the  feedback  loop  is  less  than  one.  In 
the  Nyquist  Stability  With  A  Conjecture  section,  less  conservative  conditions  for  closed-loop  stability 
are  given,  but  a  conjecture  is  needed. 

Also  in  the  PHALANX  CIWS  Aiming  Model  section,  some  design  guidelines  are  given  for  dis¬ 
turbance  rejection  and  sensitivity  to  noise.  The  disturbance  is  a  perturbation  on  the  system  output. 
Hence,  the  disturbance  is  modeled  as  an  addition  to  the  system  output.  The  noise  acts  to  corrupt  the 
measurement  of  the  output  signal.  Hence,  the  noise  is  modeled  as  an  addition  in  the  outer  feedback 
loop.  Thus,  the  filter  has  two  conflicting  tasks.  First,  it  should  block  off  the  scasor  noise.  At  the 
same  time  it  must  not  cut  off  the  feedback  loop  too  much  because  no  information  about  the  distur¬ 
bance  would  be  fed  back,  which  is  the  point  of  the  feedback  loop.  The  idea  is  to  observe  that  the  fre¬ 
quency  spectmm  of  the  disturbance  and  noise  are  different.  The  disturbance  has  a  low  frequency 
spectrum,  while  the  noise  has  a  high  frequency  spectrum.  Disturbances,  such  as  aim  bias,  would  be 
slowly  varying.  Noise,  for  instance,  due  to  thermal  effects  on  measurement  sensors,  would  have  a 
much  higher  frequency  range.  Thus,  the  desired  filter  is  a  low  pass  filter.  The  control  design  pro¬ 
cedure  used  in  this  report  is  patterned  after  the  theory  developed  in  Fruedenberg-Looze^. 

An  example  is  presented  to  illustrate  how  the  stability  conditions  and  performance  measures, 
which  are  developed  in  this  report,  are  used  to  design  a  satisfactory  control  system.  The  task  of  the 
designer  is  to  select  the  proportional  feedback  gain,  the  integral  feedback  gain  and  the  filter. 


2 


NSWCDD/TR.92/243 


PHALANX  CIWS  AIMING  MODEL 

In  this  section,  we  discuss  the  model  for  the  PHALANX  gun,  which  is  taken  from 
Bailey-Pnce  .  Figure  1  shows  the  angle  geometry  of  the  system.  Figure  2  shows  the  timing 
diagram.  In  this  diagram,  the  range  of  the  bullet  increases  with  time  while  the  range  of  the  target 
decreases  with  time.  Figure  3  shows  the  block  diagram  for  the  system,  using  Laplace  transform 
transfer  functions.  Following  is  a  list  of  variables. 

0^ ;  Aim  pointing  angle  (reference  input) 

0^ :  Measured  gun  angle 

0^ ;  Bullet  angle  at  target 

0^;  Disturbance  input 

n :  Measurement  noise  input 

tf :  Bullet  flight  time  (seconds) 

:  Estimated  bullet  flight  time  (seconds) 

This  is  a  one-dimensional  model.  In  real  life,  there  would  be  two  angular  channels,  azimuth  and 
elevation,  but  they  would  be  independent  of  each  other.  The  plant  is  modeled  by  a  pure  integrator. 
We  use  a  proportional-integral  (PI)  controller,  the  gains  are  Kp  and  Ki.  The  derivative  of  0^  is  also 
a  control  input.  This  enables  the  gun  to  track  inputs  of  a  higher  degree  than  step  inputs;  e.g.,  ramp 
inputs.  There  is  an  inner  loop  feeding  back  0^.  This  loop  takes  out  any  error  of  the  gun  aiming 
mechanism  to  the  reference  input.  The  outer  loop,  feeding  back  0/, ,  takes  out  errors  due  to  any  dis¬ 
turbances.  This  includes,  in  particular,  the  gun  aim  bias  (the  gun  aim  bias  is  modeled  as  a  distur¬ 
bance).  There  is  a  time  delay,  because  0^  is  the  angle  of  the  bullet  stream  at  the  target.  Tlie  flight 
time  of  the  bullets,  tj- ,  is  from  the  gun  to  the  target.  To  get  the  gun  aim  bias,  0^,  is  compared  to  the 
expected  bullet  stream  angle  at  the  target,  0^^  delayed  by  the  estimated  bullet  flight  time,  tf .  The  goal 
of  the  filter,  F  {s ),  is  to  filter  out  the  sen.sor  noise.  Note  that  fy-  changes  as  the  range  to  the  target 
changes. 
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FIGURE  1.  PHALANX  CIWS  ANGLE  DIAGRAM 


FIGURE  2,  PHALANX  CIWS  TIMING  DIAGRAM 
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A  point  to  consider  is  whether  this  system  is  a  time-varying  system.  At  first  glance,  ftecausc  if 
is  not  fixed,  it  may  be  fatncludcd  that  it  is  a  time-varying  system.  In  an  implcmcnlalion,  if  will 
change  continuoa.'^  as  the  utrget  maneuvers,  and  if  will  change  at  discrete  instances  as  it  is  recom¬ 
puted  by  software  within  the  gun  system.  Apparently,  there  is  no  problem  witl.  riassifyfng  the  s>'stem 
as  lime-varying.  However,  Uicre  is  another  point  of  view.  The  sy.stcm  may  also  he  elassifscd  a.s  a 
time-invariant  adaptive  .system.  If  a  target  is  engaged  one  hour  later  (or  a  day  later)  the  system  will 
perform  the  same;  hence,  it  is  time-invariant.  The  bullet  flight  time,  ij ,  ma>  be  thought  of  a.s  a  gen 
eralized  input  and  tf  is  then  an  adaptation  parameter  computed  by  system  software. 

The  first  thing  to  do  in  employing  the  Small-Gain  Theorem  is  to  compute  the  required  system 
transfer  functions.  We  need  the  transfer  functions  (9|i,/6^ )(.v )  and  )(.v  j.  These 

may  be  computed  using  Ma.son‘s  gain  ruleV  We  organize  the  loop  gaias  and  path  gaias.  There  are 
three  loops. 


L,  =- 


(Kps  +  K/) 


{KpS  Kf) 

L-,  -  - - e 


..2 


'F(.v) 


{K  s  -t-  K,) 

- £ - ^ ^  f 


There  arc  two  paths  from  9^  to  0^ . 


Pf  - 


(KpS  +  Ki) 


with  associated  loop  touching  factor  Aj^=l,  and 


with  A2  =  1.  There  is  one  path  from  9^  to  9^.  =  1  with  The  path  from  n  to 

0J,  has  gain 


-F  (s ) 


{Kps  I  Kf) 


-l,s 
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The  loop  touching  factor  is  A"  =1,  We  have  that  the  Mason  loop  gain  is  A  =  \-L  (there 

are  no  no  touching  loops).  Now, 

e7'^'  = - 5 - 


1  + 


KpS  +  Kj 


-hS 


(KpS  _  _i 

1  +  - ; - {l+f(i-)[e  ^  -e  ' 


(.5^  +  Kps  +  Kf  )(2 


-IfS 


{KpS  +  K,){  \  + 

If  tf-tf,  then  (6^/0^)(5)  becomes 

A. 

-{s)  =  e~‘^' 


1) 


(1) 


0. 


;i.e..  the  bullet  angle  at  the  target  is  equal  to  the  aiming  angle  delayed  by  seconds  (so  the  aim  bias 
is  compensated  for)  which  is  what  we  want.  Right  away,  we  see  from  Equation  (1)  that  the  time 
delay  estimation  error  is  attenuated  or  amplified  by  F.  If  0^  is  not  used  as  a  control  input,  then 
F2  =0  and  (0{,  /Og  )(s )  1  becomes 

{  KfS-HC,  1 

(S)  ^  ^ f  . 

0a  ^  -i-  KpS  -i-  Kf 

which  is  a  type  zero  system  and  so,  only  step  inputs  can  be  tracked. 

The  disturbance  transfer  function  is 


{KpS  +Kf)  ,, 

1  + - ^ - [\  -  e  '  F{s)) 

_ £2 _ 

{KpS+Kf)  : 

1+  ...  {l-hF(s)le 
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i  -  +  (Kps  +K,){  \  -  e~‘''Fis)} 
.v’  +  (Kps  +  K,){1+  Fis)\e~‘^'  - 


1  -T(5) 


(2) 


where 


Tis)  = 


The  noise  transfer  fenciion  is 


F(s){KpS  +Ki)€ 


l,S 


s^  +  iKpS  +  K,){1  +F(5)le"''" 


— /.V)  -  — _ 

n  A 


0) 


-Fis) 


(KpS  +K{) 


s‘ 


(KpS  -‘r  Kj)  , 

1  +  . r  /-{I  +F(5)fc>  '^"]) 


-F(s)(KpS  +Ki)e  ' 


t,s 


+  (KpS  +  K,)[  \  +  F{s)[e'‘''  - 

=  -7*^)  (4) 

These  transfer  functions  are  linear  with  time  delays.  Equations  (2)  and  (4)  compare  to  the  equation 
(2.3.3)  found  in  Reference  2. 

Wc  next  consider  the  stability  of  the  system.  The  first  task  is  to  assure  the  stability  of  the  inner 
loop  as  if  it  were  a  stand-alone  system.  The  transfer  function  of  the  inner  loop  is 


e. 


Kps  +  Ki 


(5) 


KpS  +  Kj 

.  or  stability,  the  poles  of  the  inner  kxtp  charactcri,stic  polynomial  +  KpS  +  !';  must  be  in  the  left 
half  plane.  This  is  done  by  selecting  Kp  >  0  and  K/  >  0.  which  also  makes  the  inner  loop  a 

e 


minimum  phase  system.  Note  that  when  or  =  yK/, 


e. 


(yto) 


>  1. 
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We  use  the  Small-Gain  Theorem  to  give  sufficient  conditions  for  the  stability  of  the  overall  sys¬ 
tem.  In  fact,  the  conditions  will  be  conservative,  which  is  a  consequence  of  the  broad  applicability  of 
the  Small  Gain  Theorem,  which  includes  time-varying  systems  with  time  delays.  Frequency  domain 
conditions  will  be  given  because  they  will  be  the  most  useful  in  conjunction  with  the  noise  sensitivity 
and  disturbance  rejection  conditions.  To  consider  the  stability  of  ihis  system,  it  will  be  convenient  to 
redraw  the  block  diagram,  Figure  4.  Because  0^  and  0^  are  inputs  to  the  system,  0^  is  restricted 
such  that  both  l0a(Ol  and  l0^(Ol  arc  bounded.  By  superposition,  the  stability  of  the  system  is 
determined  by  the  stability  of  the  system  with  switch  1  (sw#l)  closed  and  switch  2  (.sw#2)  open  and 
vice  versa.  To  use  the  Small-Gain  Theorem,  we  need  the  closed-loop  gain  for  each  of  these  cases. 
The  system  will  be  stable  if  the  closed-loop  gain  of  the  system  is  less  than  one.  (Sec  Corollary  15,  p. 
43  of  Reference  1.)  When  sw#l  is  closed  and  sw#2  is  open,  the  loop  transfer  function  is 


L  =  F{s)-{e  >  -  e  '  y 


(6) 


(5-'  -t-  KpS  -f  K,) 

When  sw#2  is  closed  and  sw#l  is  open.  Figure  4  may  be  redrawn  to  look  like  Figure  5(a).  But  Fig¬ 
ure  5(a)  is  stable  if  Figure  5(b)  is  stable.  Figure  5(b)  is  the  same  as  Figure  5(c).  These  may  be 
verified  e.g,  with  Mason’s  rule.  Hence,  the  loop  transfer  function  when  sw#2  is  closed  is  again  given 
by  Equation  (6).  By  the  Small-Gain  Thcuion,  a  sufficient  condition  for  the  stability  of  the  system  is 
that  the  loop  gain  given  by  Equation  (6)  is  less  than  one. 


The  loop  gain.  Equation  (6),  is  shown  in  Figure  6.  We  need  the  gain  from  e  to  the  point  Jt3  to 
be  less  than  one.  Bccau.se  this  is  a  linear  system,  we  just  need  to  consider  what  happens  when  a  sig¬ 
nal  with  unit  amplitude  is  input;  e.g.,  e(t)  =  cos(cOof ).  The  signal  at  point  x  j  is 
G  j  (0)(})cos(0)(/  +  (|>j)  where 


G,(0))  = 


Kpjbi)  +  Kf 
-0)^  +  Kpju)  +  Kf 


O) 


and  where  <})j  is  a  constant.  The  gain  of  e  up  to  X]  is  GjCoJo)  (i-c.,  =  Gi(tOo)).  At  ^2  the 

.signal  is  G2(ct)o)G  i(o)o)cos(a),  +()),  +  |a(/)) 
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(c) 

HGURE  5.  EQUIVALENT  BLOCK  DIAGRAMS 
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FIGURE  6.  LOOP  GAIN 


RGURE  7.  NYQUIST  PLOT  OF  (KpijOi)  +  AT/VC-CO^) 
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where 

G2i0i)  =  I =  I  ?-sin((r/  -  t^)u)/2)\  (8) 

and  where 

n(( )  =  angles =  <'/  ~  '/>“  ^g, 

is  a  time-varying  phase  shift.  Because  |X(r)  is  time  varying,  the  signal  at  X2  is  an  exponentially 
modulated  signal.  It  is  not  a  signal  with  frequency  (Oq,  rather  it  has  a  spectrum  centered  about  the 
carrier  frequency  COq.  Its  bandwidth  may  be  computed  by  Carson’s  rule^.  Let  (tj-  -  }j-)  have  max¬ 
imum  deviation,  <})^,  and  maximum  frequency,  W .  The  deviation  ratio  is  A  =  ®y  Carson’s 

rule,  the  bandwidth  of  the  signal  at  X2  is  5(cOo)  =  2  (A  -t-  O  W.  We  will  assume  that  the  spectrum 
of  X2  is  entirely  contained  within  the  bandwidth  S(O)o)  on  either  side  of  the  carrier  frequency  coo. 
The  gain  of  the  signal  at  X2  is  G2(cOo)  G  [((Hq),  the  same  as  it  would  be  if  }i.(r)  were  a  constant.  We 
compute  the  gain  atX3  due  to  e{t)  =  cos(tOo/)  as  follows: 


Compute  B  (cOq)  and  find 


G3(5(o)o))=  max  lF(yo))l 

(Ho  -  B  (tOo)  <  CO  1  oJo  +  B  (cod) 


then  the  gain  at  x  3  is  bounded  by 


G (coq)  -  G-^iB  (o)o))  G2(0)o)  G  i(coo)  (l  l) 

The  loop  gain  for  the  overall  closed  loop  system  is  bounded  by 

G  =  sup  G((Do)  fj2) 

coo 

We  take  a  brief  look  at  the  size  of  the  loop  gain.  For  higher  frequencies,  looking  at  Equation 
(11),  things  should  be  fine  because  Gi(to)  goes  to  zero,  G2(to)  is  less  than  two  and  the  filter  is 
selected  such  that  G3(a))  is  bounded.  At  low  frequencies,  G  i(a))  is  about  one  and  G2((o)  approaches 
zero.  The  product  of  G](to)  and  Gj(B  (to))  must  be  less  that  one  half  before  G2((0)  reaches  two  for 
the  first  time.  For  faster  inner  loop  dynamics,  Gj(co)  will  break  at  a  higher  frequency,  so  that 
I  ty  -  I  must  be  .smaller;  hence,  a  more  accurate  is  required.  Suppose  that  the  filter  is  a  low 
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pass  filter  of  the  form  F(,s)  =  A/(s  +  A ).  Again,  looking  at  Equation  (8),  if  the  filter  breaks  at  a 
higher  frequency,  tf  must  be  more  accurate.  Also,  for  this  type  of  filter,  if  W  is  larger  G  3(6  (to)) 
will  drop  at  a  higher  frequency,  so  that  a  more  accurate  tf  is  required.  (See  Equation  (10).) 

Consider  the  disturbance  and  noise.  As  mentioned  in  the  Introduction,  disturbances  are  a  factor 
at  low  frequencies,  while  noise  is  a  factor  at  high  frequencies.  Looking  at  Equations  (2)  and  (4),  dis¬ 
turbance  rejection  (at  low  frequencies)  and  low  noise  sensitivity  (at  high  frequencies)  are  achieved  if 
II  r(/o))ll  is  small  at  high  frequencies  and  T(s)  =  1  at  low  frequencies.  If  the  system  is  run  open- 
loop:  i.e..  F  (5 )  =  0,  the  system  will  not  be  susceptible  to  noise;  however,  we  need  the  feedback  to 
reduce  the  impact  of  the  disturbance  to  the  system.  At  low  frequencies  ~  1  and  because 

-  tf)  ~  0,  we  have  that  -  e  =  0.  Hence,  at  low  frequencies 


r(jcD)  = 


FO(o)(A^/>yto-b  Ki) 
—03^  -f  {K-p _/o)  +  Kf) 


=  FU(^) 


so,  F(/C0)  close  to  one  at  low  frequencies  gives  better  disturbance  rejection.  Agai.n,  suppose  the  filter 
is  of  the  form  Fis)  A/(s  +  A).  We  see  from  Equations  (2)  and  (4)  that  if  F  breaks  at  a  higher 
frequency  there  is  better  disturbance  rejection,  while  F  breaking  at  a  lower  frequency  reduces  sensi¬ 
tivity  to  noise.  A  design  example  will  be  presented  later  to  clarify  this  discussion. 


NYQUIST  STABILITY  WITH  A  CONJECTURE 

The  Nyquist  stability  conditions  will  be  used  in  this  section  to  analyze  the  closed-loop  stability 
of  the  system.  To  handle  the  time-varying  nature  of  this  syslctn,  the  following  conjecture  will  oe 
accepted.  Let  T  be  large  enough  so  that  all  conceivable  tf,tf  are  in  (0,7].  If  the  system  is  stable 
for  all  fixed  (time  invariant)  tf.'tf  e  [0,7],  then  the  system  is  stable  for  lime  varying  tf,  'tf  so  long 
as  they  remain  in  the  interval  [0,7],  We  do  not  mean  that  this  conjecture  could  be  proven  (as 
counter-examples  exist),  rather  we  want  to  find  some  conditions  which  would  be  necessary. 
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Equation  (1)  has 


e. 


-ffS 


<s) 


(s^  +  KpS  +  K,)e 


= 


s^  +  (KpS  +^/){l  +F(s)[e~‘'^ 

(KpS  +  Ki)F(s)[e'^‘'*‘^'^ 

+  (KpS  +  K,){  \  +  Fis)[e~‘'^  - 

{KpS  +K,)  ^  ^  _2;,  ,  j 


+  F(s) 


(KpS  +  K,)  _} 

1  + -  - {1+F(s)[e  '  -  e  '  ]} 


(13) 


If  tj-  =  /y,  the  second  term  in  Equation  (13)  would  be  zero  and  (0^/0g)(5)  would  be  a  stable  transfer 
function.  Consider  the  stability  of  the  second  tenn.  We  need  the  zeros  of 


(Kp  s  +  Kf)  - 

1  +  -  -{1  +Fis)(e~‘f'  -e  ''")) 


(14) 


to  be  in  the  left  half  plane.  First,  the  Nyquist  plot  of  the  first  factor  (KpS  +  Ki)Js^  is  drawn.  See 
Figure  7.  The  phase  margin  is  ^  and  there  is  infinite  gain  margin.  The  second  factor  of  Equation 
(14)  is  of  the  form  (1  +  e)  with  e  =  F  (s  )(e  -  e  *■'*).  If  this  term  has  a  phase  angle  of  less  than 

<|>,  the  closed-loop  system  will  be  stable.  A  sufficient  condition  for  this  is 


I  sin  *(e)l  <  ^ 


(15) 


Equation  (15)  presents  a  tradeoff  between  the  error  in  the  flight  time  estimation  error  and  the  filter 
F  (s  ).  That  is,  wider  filter  bandwidth  requires  more  accurate  tr . 


DESIGN  EXAMPLE 

The  results  presented  in  the  PHALANX  CIWS  Aiming  model  section  are  illushated  by  an  exam¬ 
ple.  The  inner  loop  stability  and  performance,  overall  system  stability,  disturbance  rejection  and  sensi- 
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dvity  to  noise  arc  considered.  First,  proportional  and  integral  gains  Kp  and  K[  are  selected  for  inner 
loop  stability  and  performance.  We  select  Kp  and  /f/  so  that  the  two  inner  loop  closed-loop  poles  are 
evenly  spaced  on  a  half  circle  of  radius  10  in  the  left  half  plane.  This  gives  Kp  =  17.33  rad  and 
Kj  =  100.0  rad/sec.  These  poles  determine  the  reaction  performance  of  0^  to  0^ . 


We  choose  the  Alter  to  be  a  Arst-order  low  pass  Alter, 


F(5)  = 


A 

is  +A) 


(16) 


The  pole  of  F  (j )  is  selected  so  that  the  Alter  breaks  at  a  higher  frequency  than  the  disturbance  fre¬ 
quency  band  and  at  a  lower  frequency  than  the  noise  frequency  band.  With  this  choice  of  F(s), 
Equation  (3)  becomes 


....  .  AiKfS +Kiye-^'’ 

T(s)  = - = — 

(5^  KpS  -I-  K,yis  +  A)  +  A  (KpS  +  Ki)  [e~‘f'  - 


(17) 


and  Equation  (6)  becomes 

A-ie~‘'^  -e'‘'^)iKpS  +  K.) 

I  = - ^ .  (Ig) 

(5  ^A)is^^KpS  ^-K,) 

We  want  to  compute  the  maximum  possible  bullet  iight  time  estimation  error  that  maintains  sta¬ 
bility  of  the  overall  closed-loop  system.  Note  that  (ty-  -  ty)  enters  into  Equation  (11)  in  two  places. 
The  difference  (/y  -  t/)  appears  in  G2((Oo)  (See  Equation  (8).)  and  in  GjiBiOg))  via  the  bandwidth 

B  (cOq)-  (Thus  B  (to)  =  2  (-^maxl  ty  -  ty  I  o)  -f  I)  W.)  We  set  values  for  A  and  say  A  =1.0 


rad/sec  and  W  =  1.0  rad/sec.  With  this  value  for  W ,  we  And  the  maximum  I  ty  -  ty  I  such  that  the 
loop  gain  is  less  than  one.  Figure  8  shows  loop  gain  vs.  frequency  for  I  /y  -  ?y  I  =  0.29  sec.  The 
gain  is  less  than  one  for  all  frequencies;  hence,  by  the  Small-Gain  Theorem,  the  overall  closed  loop 
system  is  stable  (for  I  ty  -  ?y  I  <  0.29). 


As  stated  earlier,  if  A  is  larger  (the  Alter  breaks  at  a  higher  frequency),  fy  must  be  more  accu¬ 
rate;  while  if  A  is  smaller  ty  may  be  less  accurate.  Also,  a  larger  W  will  require  a  more  accurate  ty ; 
while  a  smaUer  W  allows  a  less  accurate  ty .  Several  test  cases  arc  given  in  Table  1. 


16 


I  (  )(;■«)  •  LOOP  GAIN 


NSWCDD/TR-92/243 


10"^  ^ - ‘ - ' ‘ — — ' — ' — - ‘ - - ' ^ — ‘ — 1  J - ■ — ‘ .j i 1„ 

10-2  10-1  IQO  IQl  102 

FREQUENCY  (RAD/SEC) 

FIGURE  8.  LOOP  GAIN  VS.  FREQUENCY 


FREQUENCY  (RAD/SEC) 

FIGURE  9.  DISTURBANCE  GAIN  VS.  FREQUENCY  FOR  tf  =  5.0  SEC,  r  y  =  5.25  SEC 
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TABLE  1.  MAXIMUM  ALLOWABLE 
\tf  -tf\  FOR  STABILITY 


Maximum  1  fy  ~  fy  1 

A  rad/sec 

1.0 

0.1 

1.09 

10.0 

0.1 

0.129 

1.0 

1.0 

0.29 

10.0 

1.0 

0.10 

100.0 

1.0 

0.069 

0.1 

10.0 

0.32 

1.0 

10.0 

0.068 

10.0 

10.0 

0.068 

Finally,  we  take  a  look  at  the  disturbance  and  noise  transfer  function  frequency  plots.  For  the 
remainder  of  this  report  A  =1.0  rad/sec  and  W=1.0  rad/sec.  The  disturbance  and  noise  transfer  func¬ 
tions,  Equations  (2)  and  (4),  respectively,  will  be  treated  as  though  tj-  and  tj-  are  constants.  Time 
varying  ty  and  will  widen  out  the  disturbance  and  noise  spectrums.  Figure  9  shows 
)  (06/Bj/)C/O))I  vs.  0)  for  /y  =  5.0  sec  and  ?y  =  5.25  sec.  This  is  about  the  longest  flight  time  of 
the  bullets.  For  frequencies  less  than  0.01  rad/sec,  this  transfer  funcUon  has  gain  less  than  0.1; 
hence,  the  low  frequency  disturbances  are  filtered  out  and  have  little  effect  on  the  output.  Figure  10 
show.s  I  (0^ /«)(/■  (0)1  vs.  0)  for  the  same  /y  and  /y.  For  frequencies  greater  than  40  rad/sec,  the 
gain  of  this  transfer  function  is  less  than  0.01;  hence,  the  high  frequency  noise  is  filtered  out.  Finally, 
we  check  the  disturbance  and  noise  transfer  functions  for  a  short  bullet  flight  time.  Figures  11  and  12 
show  I  {0j,/0^)(/(O)l  and  I  (0^ /rt  )(y  (O)  I  vs.  frequency,  respectively,  for  /y  =  1.25  sec  and 
(y  =1.0  sec.  These  represent  approximately  the  shortest  bullet  flight  time.  Again,  the  low  frequency 
disturbances  and  high  frequency  noise  are  filtered  out.  Of  course.  Table  1  may  be  used  the  other  way 
around;  i.e.,  given  a  maximum  I  (y  -  (y  I ,  A  and  VV'  may  be  selected. 

Looking  at  Figures  9  and  11,  it  seems  that  1  (0^/0<^)(/O^)l  drops  below  0.1  at  too  low  of  a  fre¬ 
quency,  We  would  like  to  get  the  curves  in  these  figures  to  drop  down  at  higher  frequencies.  There 
is  a  physical  limitation,  however.  From  Figure  3,  it  is  seen  that  a  variation  in  the  disturbance  cannot 
be  immediately  detected  -  there  is  a  time  delay  of  ry  seconds.  This  can  also  be  seen  from  Equations 
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10-2 


]0'2 - ^ - ' ' - , ; 1 1 — - x_  . . 

10-2  10-1  IQO  10' 

FREQUENCY  (RAD/SEC) 

FIGURE  10.  NOISE  GAIN  VS.  FREQUENCY  FOR  tf  -  5.0  SEC.  i  /  = 


HGURE  11.  DISTURBANCE  GAIN  VS.  FREQUENCY  FOR  tf  =  1.25  SEC, 


102 

5.25  SEC 


UP 

}  f  =  1.0  SEC 
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(2)  and  (17).  For  I  (0j,/0^)(/'Co)l  to  approach  zero,  T(j  O))  must  approach  one  in  magnitude  and 
angle.  The  factor  in  Equation  (17)  adds  phase  lag.  We  further  illustrate  this  situation  in  fig¬ 

ure  !3.  There  arc  three  curves,  the  first  is  il  ~  (KpiJ  o))  +  Kj  )  A : 
((-CO^  +  KpiJOi)  +  Ki)'(J(33  +  /4  ))l  vs.  0).  This  is  a  plot  of  the  disturbance  gain  with  tf  =  ij  =  0 
(and  the  other  parameters  as  before.)  This  plot  represents  the  ideal  case  of  disturbance  rejection,  from 
the  point  of  view  of  the  physical  limitation;  i.e.,  the  aim  bias  is  in.stanianeously  fed  back.  If  Kp  and 
A  were  increased,  the  break  would  be  at  a  higher  frequency.  The  second  curve  in  Figure  13  is 
1  1  -  vs.  CO  with  ij-  =  2.5  sec.  This  is  an  average  bullet  flight  time.  This  represents  the 

physical  limitation  due  to  time  delay  of  how  fast  the  disturbance  transfer  function  may  approach  zero 
as  CO— >0.  The  third  curve  is  a  plot  of  I  (0j,/0^)(7to)l  vs.  co  with  ij-  -Ij-  =  2.5  sec.  so  that  there  is 
no  estimation  error.  This  curve  has  the  combined  effects  of  the  first  two  curves.  Hence,  the  third 
curve  breaks  at  a  lower  frequency  than  either  of  the  first  two  curves.  The  first  two  curves  may  be 
thought  of  as  representing  the  limit  of  the  disturbance  rejection  performance  of  the  gun  aiming  control 
system. 
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CONCLUSION 

The  analysis  presented  in  this  report  provides  new  insights  into  the  stability  characteristics  of  the 
aim  bias  compensation  function  of  the  PHALANX  CIWS  Block  1.  Also,  tradeoffs  between  the  com¬ 
peting  design  goals  of  stability,  disturbance  rejection  (e.g.  aim  bias  compensation)  and  sensitivity  to 
feedback  noise  are  presented. 

The  aiming  control  system  is  a  linear  system  with  a  time-varying  time  delay.  The  aiming  control 
system  features  an  inner  feedback  loop  servomechanism  to  aim  the  gun  and  an  outer  loop  to  null  out 
any  aim  bias.  The  time  delay  results  from  feeding  back  the  bullet  angle  at  the  target  through  the  outer 
loop.  Since  the  bullet  flight  time  changes,  the  time  delay  is  time-varying.  The  Small-Gain  theorem 
may  be  ased  to  give  sufficient  stability  conditions.  The  Small-Gain  theorem  and  the  noise  and  distur¬ 
bance  transfer  functions  are  used  to  analizc  design  tradeoffs.  Faster  inner  loop  dynamics,  a  noise  filler 
with  wider  bandwidth  or  a  higher  frequency  in  the  time  delay  estimation  error  signal  will  result  in  a 
less  stable  system.  The  filter  in  ilie  outer  loop  regulates  the  disturbance  rejection  and  noise  filtering 
bandwidlhs.  Since  the  noise  and  disturbance  transfer  function  must  sum  to  one,  more  noise  filtering 
results  in  less  disturbance  rejection  and  vice  versa.  Finally,  we  have  found  that  the  time  delay 
represents  a  limit  to  the  disturbance  rejection  bandwidth  of  the  control  system. 

The  insights  and  design  tradeoffs  can  be  used  to  screen  alternatives  in  the  first  stages  of  concept 
development  cf  new  applications  or  upgraded  versions  of  the  PHALANX  system.  A  judicious  combi¬ 
nation  of  analytical  and  experimental  analysis  will  allow  a  design  to  be  achieved  with  good  balance 
between  .stability  and  responsiveness  with  a  reduced  amount  of  lime  and  effort. 
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